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The model of cylindrical gravitational waves is employed to work out and check a recent proposal in 
Ref. [TH how a diffeomorphism-invariant Hamiltonian dynamics is to be constructed. The starting 
point is the action by Ashtekar and Pierri because it contains the boundary term that makes it 
differentiable for non-trivial variations at infinity. With the help of parametrization at infinity, 
the notion of gauge transformation is clearly separated from that of asymptotic symmetry. The 
symplectic geometry of asymptotic symmetries and asymptotic time is described and the role of 
the asymptotic structures in defining a zero-motion frame for the Hamiltonian dynamics of Dirac 
observables is explained. Complete sets of Dirac observables associated with the asymptotic fields 
are found and the action of the asymptotic symmetries on them is calculated. The construction of 
the corresponding quantum theory is sketched: the Fock space, operators of asymptotic fields, the 
Hamiltonian and the scattering matrix are determined. 
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I. INTRODUCTION 

A deep problem in quantum gravity is the dependence 
of the constructed quantum theory on the choice of co- 
ordinates. A famous example of a quantization method 
based on such a choice of gauge is the Arnowitt, Dcser 
and Misner reduction Q . There have been ideas of how 
this problem could be avoided, such as Dirac's opera- 
tor constraint equations Bergmann's work on Dirac 
observables 0], BRST method Q, Euchdean quantum 
gravity etc. None of these methods has as yet been 
really successful because new problems always emerged. 

In this paper, we focus on the method based on the 
Dirac observables. In a sense, this is the most straightfor- 
ward one: all variables that are used have to be gauge in- 
variant. However, changes of coordinates include changes 
of time so that the gauge invariance entails the time in- 
dependence: Dirac observables must be integrals of mo- 
tion. Two problems that are related to this have been no- 
ticed already by Bergmann Q: frozen dynamics" and 
scarcity: the dynamics of Dirac observables is trivial — 
they just stay constant, and it is difficult to find any 
such quantity in general relativity. More recently, the 
non-locality of such variables has been proved 0: the ex- 
pression for any Dirac observable in terms of local fields 
must contain derivatives of all orders. 

As for the "frozen dynamics" , Rovelli's idea of "evolv- 
ing constants of motion" 7] has shown one way out of the 
problem. Another, but not completely unrelated way, is 
based on the observation that any Hamiltonian formula- 
tion of dynamics needs a frame [8| . A possible frame for 



a non-trivial Hamiltonian dynamics of Dirac observables 
has been specified in and ^9 under the assumption 
that there is a symmetry. In [ij, it is shown that even 
an asymptotic symmetry suffices, at least for a simple 
finite-dimensional model. 

Concerning the scarcity, there have been different pro- 
posals of how Dirac observables could be constructed 
(e.g., 0). Also, in the case of asymptotically ffat mod- 
els, there are even complete sets of natural constants of 
motion, namely the in- and out-fields. These quantities 
have been described by DeWitt within his covariant 
perturbation theory and by Ashtekar in general. 

Finally, the non-locality of the kind found need not 
always lead to a real problem. For example, in any non- 
trivial quantum field theory in Minkowski spacetime, the 
expressions for the in- and out-fields in terms of the local 
fields at a finite time are also badly non-local. However, 
such expressions are not needed. What is needed are 
the expressions of the in- and out-fields in the in- and 
out-regions, respectively, which are local. 

In the present paper, we are going to strengthen the 
point made in favour of the Hamiltonian frame based 
on asymptotic symmetries. We study a field model, ex- 
tending thus the cases in which the idea works to infinite- 
dimensional systems (see also [l^ ) . The simplified model 
that we investigate consists of cylindrical gravitational 
waves, also called sometimes Einstein-Rosen waves (see, 
e.g., and ^3)- 1^ seems that this model can be 
used as an example of almost everything. Thus, Kuchaf 
has studied the embedding variables in canonical theory 
employing the waves |l8j |. Torre has expressed a com- 
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plete set of Dirac observables in terms of local fields [l9| . 
Ashtekar and Pierri have found a considerable simpli- 
fication of Hamiltonian dynamics of the waves using a 
suitable gauge |20l |. 

Most important, in pol | an asymptotic boundary term 
analogous to the ADM energy has been added to the ac- 
tion for the first time. This is an achievement because 
the cylindrical wave spacetimes are not asymptotically 
flat due to the cylindrical symmetry so that the well- 
known results in general relativity cannot be used. The 
model can, nevertheless, be reduced from four to three 
dimensions by removing the direction of the translational 
part of the cylindrical symmetry. The resulting system 
has the form of gravity coupled to a scalar field and the 
geometry of the three-dimensional spacetime is asymp- 
totically flat in certain sense (locally asymptotically flat). 
This has been shown in and [22| , where the boundary 
term has been found. Our analysis is, therefore, based 
on the results of Ashtekar and Pierri. 

In Sec. 2, the model of free Newtonian particle moving 
in one-dimensional space is used to explain how a sym- 
metry provides a frame for a non-trivial Hamiltonian dy- 
namics of integrals of motion. The relevant notions and 
relations are introduced. The method of reduction by a 
choice of gauge is also described within this framework. 

Sec. 3 summarizes the well-known results on the cylin- 
drical waves that we shall need, focusing on the asymp- 
totic properties. Sec. 4 reviews and modifies the results 
of Ashtekar and Pierri so that they become compatible 
with the theory of Sec. 2. The meaning of asymptotic 
symmetry and asymptotic time are explained and their 
relation to "ordinary" symmetry and time is clarified. 
The way in which the asymptotic symmetry and time 
determine a frame of a Hamiltonian dynamics of Dirac 
observables for the cylindrical waves is described. An 
important new point made in Sec. 4 is a clean sepa- 
ration between gauge transformations and symmetries. 
In |23|, [l^l and 'll], the group of diffeomorphisms has 
been divided into gauges and symmetries according to 
the action of its elements at infinity. However, it has 
never been completely clear where the boundary is to be 
drawn. When we tried to apply this idea to the cylin- 
drical waves, some strange paradoxes have appeared. It 
has turned out that the action must be parametrized at 
infinity similarly as in the case of geometrodynamics of 
Schwarzschild black holes privileged spacetime co- 
ordinates at infinity must be added to the set of canonical 
coordinates of the phase space. Then, all gauge transfor- 
mations including reparametrizations at infinity are gen- 
erated by constraints, while all asymptotic symmetries 
are generated by expressions in the momenta conjugate 
to the privileged coordinates at infinity. 

The Poisson algebra of the asymptotic Dirac observ- 
ables listed in Sec. 3 is calculated in Sec. 5. The physical 
phase space is defined by a complete set of Dirac observ- 
ables and their Poisson brackets. In Sec. 6, the action of 
the asymptotic symmetries on the physical phase space 
is found and the canonical generator of the continuous 



group of asymptotic time translations is written down. 
The corresponding Hamiltonian dynamics in the physi- 
cal phase space is described. It may seem paradoxical 
that all these gauge-invariant structures (complete set of 
Dirac observables, their Poisson brackets and the action 
of asymptotic symmetry on them) can be and, indeed, 
have been calculated from the gauge-dependent action 
obtained by Ashtekar and Pierri after a choice of gauge 
(cf. also [13 )■ The justification thereof is given in Sec. 4. 

Finally, Sec. 7 gives a brief sketch of how the results 
of Sees. 5 and 6 can be employed for one of possible con- 
structions of quantum theory. The Fock space and the 
basic operators on it are specified. The Hamiltonian op- 
erator and the scattering matrix are determined. 



II. DYNAMICS OF DIRAC OBSERVABLES 
IN A ONE-DIMENSIONAL MODEL 

A finite-dimensional example can illustrate the main 
aspects of our method. The underlying geometric frame- 
work has been developed for general finite-dimensional 
systems in Refs. [H, 0, |2| and [lOl. 

Consider a free Newtonian particle of unit mass in one- 
dimensional space evolving in the Newtonian time T. Let 
its position be denoted by Q. The system is described by 
the action 



s[Q] = \ dT (g 



Ti 

This action is invariant under the translation 



T 



(2.1) 



(2.2) 



for any time T and any real parameter r. This is a one- 
dimensional group of symmetry that will play a crucial 
role in what follows. The symmetry implies via Noether 
theorem the conservation of energy 



E 



1 



The corresponding Hamiltonian action is 



S[Q,Pq] = I ' dT(pQQ 



1 . 



(2.3) 



The one-form Q — Pq dQ contained in the action is called 
the Liouville form. The space with the coordinates Q and 
Pq is the physical phase space Fi ; it carries the symplectic 
two- form ill = d6 = dPQ A dQ. 

The equations of motion implied by action H2.3|l are 



Q.T - Pq 
Pq,t = . 

Their general solution is 



Q{T) = q- 

Pq{t) = p 



pT , 



(2.4) 
(2.5) 



(2.6) 
(2.7) 
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where {q,p) G are constant for a particular solution. 
We can, therefore, define the space of solutions T2 as 
with coordinates q and p. An important observation 
is that r2 also carries a symplectic structure. Indeed, 
Eqs. H2.6() and (|2.7|l can be considered as a T-dependent 
canonical transformation because they imply 



PgdQ ^ pdq + d 



{Q-q? 
2T 



where, of course, the action of the differential d "hits" 
only the variables Q and q, not T; {Q — g)^/2T is the 
generating function. The resulting symplectic form Q2 
of r2 is dp f\ dq. The dynamical change 5d within the 
time increment 5T in r2 is trivial: 



(5dP = , 5dq = 



(2.8) 



We can also consider the solutions H2.6|l and (|2.7|) as 
defining two scalar fields on the one-dimensional time 
manifold R with coordinate T, the so-called background 
manifold. The push forward by the symmetry transfor- 
mation 1)2.2(1 acts on these fields as follows 



Pq{t) 



Q{T- 

Pq{t- 



r) , 



(2.9) 
(2.10) 



This symmetry action preserves, of course, the property 
of being a solution of the equations of motion; if q and 
p describe the untransformed solution and q' and p' the 
transformed one, then Eqs. H2.6|l and (|2.7|l yield 



q = q-pr 



p ^p . 



(2.11) 



This action of the symmetry transformation on the space 
r2 is canonical and is generated by the function —p^jl 
(as the Poisson brackets show): 

hq - -p<5T = {g,-yp2| ^ (2.12) 

bsP - = {p,-yp2| ^ (2.13) 

where 6^ is the change caused by the symmetry transfor- 
mation. 

The key observation is now the following. If we com- 
pare the time change in q and p due to the dynamics (Eq. 
(12.81) 1 with that due to the symmetry ((Eqs. H2.12|l and 
(12.13(1 1. we find that the relative time changes Sdq — Ssq 
and SdP — SsP are formally identical to the original dy- 
namics generated by the Hamiltonian ^q/2 in Fi. 

The description of the dynamics of our model can be 
made generally covariant by parametrization (see, e.g., 
[22,]). An arbitrary time t = t{T) is introduced and the 
action ((2.1(1 is brought into the constraint- Hamiltonian 
form (i.e., the Hamiltonian is a linear combination of 
constraint functions — a typical property of generally co- 
variant systems, see 0): 



S[Q,Pq,T,Pt;N] 



(2.14) 



dt[PQ Q,t + PrTt ~N{Pt + ^ P|) 



Let the space P consist of all initial data Q, Pq, T and 
Pt for the field equations implied by action ((2.14(1 . The 
space V equipped with the symplectic form fl-p derived 
from the Liouville form of action (12.14(1 . 

n-p = dPg /\dQ + dPr A dT , 

is called the extended phase space of the system. is a 
Lagrange multiplier and C = is a constraint with 



C^Pt 



1 



P 



Q 



(2.15) 



being the constraint function. The surface C inV defined 
by the constraint is called the constraint surface. Since 
equation C — can be solved for Pt, we can choose the 
functions Q, Pq and T as coordinates on C. 

The transformations that are canonically generated 
by the function H{N) = NC can be considered as 
reparametrizations — general changes of the time param- 
eter. They are, therefore, analogous to gauge transfor- 
mations. The corresponding gauge group acts along the 
constraint surface C. At the same time, its orbits in C 
coincide with the solutions ((2.61) and (12.7(1 . The physical 
phase space r2 of gauge-equivalent solutions can, there- 
fore, be identified with the quotient of the constraint sur- 
face by the gauge orbits: 



C 
OA 



(2.16) 



In coordinates {q,p) on r2 and {Q, Pq,T) on C, the pro- 
jection Proj^g^p^-) derived from ((2.16() is 

[q,Pq,t) eC ^ (g,p) = {Q-PqT,Pq) eT2 . 

(2.17) 

The symplectic form D,2 on r2 can be obtained from fl-p 
as follows. First, fip is pulled back from V to C. This 
yields a two-form fie degenerated along the gauge orbits 
in C. The projection Proj((^ _^ p^) ((2.17(1 determines the 
symplectic form 172 as the unique solution of the equa- 
tion Qc ~ -Pi'oj(c^r2) ^2 where * denotes the pull-back 
mapping. 

We shall also need the concept of transversal surface 



T dC. 

that 



This is a section 



Idr 



C in the sense 



(2.18) 



with respect to the projection ((2.17() . Whenever C ad- 
mits such a section, the surface T — (7(r2^c)(r2) is a 
copy of the physical phase space r2. A bijection between 
any surface T and the physical phase space F2 can be 
defined by restricting the projection Projj^^p^) to each 
particular T. The symplectic form fl2 induces through 
each such bijection a unique symplectic form ^x- In this 
way each T also becomes a phase space with a symplectic 
structure that is isomorphic to that of F2. Transversal 
surface T is called regular if if it is not tangential to gauge 
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orbits at any of its points. If T is defined by the equa- 
tion F{Q, Pq,T) — 0, then the regularity condition is a 
non- vanishing Poison bracket 

{F,niN)}\c^O yN . (2.19) 

The meaning of the regularity of transversal surfaces sim- 
ply is that the gauge condition breaks the gauge com- 
pletely. Systems that are not pathological possess many 
transversal surfaces. 

One way to quantize a generally covariant system is to 
reduce it to an unconstrained system of the kind H2.3|l . 
We shall now show two methods of reduction: by a gauge 
condition and via Dirac observables and symmetry. 

A. Reduction by a choice of gauge 

A gauge condition is a choice of a particular family of 
regular transversal surfaces that foliate C. Let the family 
be given by the set of equations 

F{Q,PQ,T)^f- (2.20) 

for each fixed real T, one surface Tf of the family is 
defined. The reduction using the condition H2.20|l can 
proceed as follows. Suppose that a canonical trans- 
formation in V is known between the original vari- 
ables {Q,Pq,T,Pt) and some canonical coordinates 
{q, p, T, P) that have been chosen so as to contain T. 
Then, using the regularity condition of the gauge, one 
can show that the constraint C = is solvable with re- 
spect to P and so can equivalently be written as 

P + 7Y(f,g,p) = 0, 

where Ti. is some smooth function. The canonical trans- 
formation brings action (|2.3|) to the form 

S{q, p, f,P,N)^ r dt [Pit + FT t - N{P + H)] , 
Jti 

where iV is a new Lagrange multiplier defined by 

NC = N{P + ii) . 

Next, T is chosen as the integration variable t and the 
action is restricted to the constraint surface. The result 
is 

S{q,p)^ df{pq f^n) . (2.21) 

By this, the reduction is finished. 

A problem with this kind of reduction is that the new 
variables g, p as well as the new time T are not the same 
as the original variables Q, Pq and T. Classically, the 
two actions (I2.21|) and are equivalent, because they 
are related by an extended gauge transformation. The 



two quantum mechanics, however, that are obtained by 
the standard quantization method from them, cannot be 
unitarily equivalent [23|: the transformation H2.2Uf) be- 
tween the respective times involves operators, while each 
of the times must be a parameter in the respective quan- 
tum mechanics. 



B. Reduction using Dirac observables and 
symmetry 

A Dirac observable o{V) is a function o : V whose 
Poisson bracket with the constraint C vanishes when re- 
stricted to C. Dirac observables are gauge invariant. 

The correspondence between Dirac observables on V 
and functions on r2 is the following: Each function / : 
r2 i-^ M determines a function /o Proj^^^p^-j on the 
constraint surface via the projection mapping H2.17|l . In 
the chart (Q, T, Pq) on C such a function has the form 

fiQ^PQT,PQ). (2.22) 

The next step is to extend this function from C toV. Let 
us denote such an extension by o : P i— > R. The only 
condition is that o be a smooth function on V the values 
of which coincide with f{Q — PqT, Pq) at C. In this way, 
the original function / on T2 defines an equivalence class 
{0} of functions a onV. One can show that any two ele- 
ments oi and 02 of {0} satisfy 02=01+ NC, where N is 
a smooth function on V (see, e.g., 0). It also follows im- 
mediately from the construction that o is constant along 
orbits so that it is a Dirac observable. Thus, a whole 
class of Dirac observables corresponds to one function on 
r2 (one often speaks about Dirac observables meaning 
these classes). 

On the other hand, each Dirac observable o defines via 
the restriction to C a function on C that is constant along 
orbits. Such a function determines, in turn, a unique 
function / on r2. 

The Poisson brackets between Dirac observables can 
be calculated using the symplectic structure of the ex- 
tended phase space V. It is easy to show that the 
Poisson bracket {01,02} of two Dirac observables is again 
a Dirac observable and that the Poisson brackets {oi, 02} 
and {oi+NiC, 02 + N2C} lie in the same class. Thus, the 
Poisson brackets between the equivalence classes {0} are 
well defined. Moreover, if the class with the representa- 
tive Oi corresponds to the function /i, i = 1, 2, on r2, and 
the class with the representative {01,02} corresponds to 
/, then 

{fi, /2}r2 = / , 

as is shown in Ref. 9]. It follows from this that a com- 
plete set of Dirac observables, together with their Poisson 
algebra, determine the structure of the physical phase 
space r2 . A complete set of Dirac observables separates 
gauge orbits, and in simple cases can be used as a coor- 
dinate system on the quotient space C/Orb = r2. For 
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our model, a complete set is formed by the functions 
oo = Q - PqT + NqC and oi ^ Pq + NiC, where No 
and iVi are smooth functions on P. A simple calculation 
gives that the only non-trivial bracket is 

{00,01} -l + A^C, 

where N — {Nq, 01} + {oq, Ni}; the Dirac observables oq 
and oi correspond to the functions q and p on 

The symmetry group (|2.2|l acts on the extended phase 
space as follows 

iQ,PQ,T,PT)^iQ,PQ,T + T,PT) , (2.23) 

and is, therefore, generated by the momentum Pt conju- 
gate to T. Observe that Pt itself is a Dirac observable; 
one can prove that any continuous symmetry group 
is generated by a Dirac observable. Now, Pt has a non- 
trivial action on Dirac observables. For example, 

{00, Pt} = -oi+NoC , 
{01, Pt} = NiC , 

where Nq and Ni are suitable functions on V. The change 
of Dirac observables referred to the symmetry as "zero 
motion" is, therefore, non-trivial. It is easy to see that 
this change is generated by the function —Pt- The value 
of Pt at C is, however, 

-Pt\c = \PI , 

and it lies in the class o^/2 + NC . The corresponding 
function on r2 is, therefore, H = p^/2, and it plays the 
role of the Hamiltonian of the constructed dynamics. In 
this way, we have recovered the dynamics and the phase 
space Fi of the original system so that the reduction is 
accomplished. 

Mathematically, any symmetry of V that is not pure 
gauge transformation can generate a non-trivial evolution 
of Dirac observables on F2 because it defines a non-trivial 
mapping between gauge orbits in C. By projection to F2 
a symmetry is obtained which can be interpreted as the 
generator of a dynamical evolution on F2. Physically, 
it must be additionally required that the symmetry be 
privileged by the situation at hand. Only then its role 
as a true Hamiltonian on F2 can be justified. Here, the 
constant translation (|2.23|) is physically privileged by the 
arguments leading from H2.1() to (|2.14l) and in particular 
by the fact that it yields through Noether's theorem the 
energy of the Newtonian particle in the privileged refer- 
ence system T. The transformation H2.23|l is indeed a 
symmetry of V because the Poisson bracket of its gener- 
ator, — Pt, with the constraint function l|2.15|l vanishes 
on C. 

The following observation is very important. If Pt 
generates a symmetry that leads to the Hamiltonian H 
in F2, then so does Pt + N'C for any smooth N': the 
dynamics of Dirac observables is uniquely determined by 
the whole class of symmetry generators. Why is this 



important: In our simple model, we have a unique sym- 
metry and it is generated by Pt- The reason is that our 
model is a so-called "already parametrized system" with 
a privileged time T . Indeed, there also is a privileged 
choice of gauge due to this fact: P(Q, Pq,T) = T, which 
leads to the "right" action H2.3|) by the reduction proce- 
dure of Sec. 2.1. However, many models of real interest, 
such as general relativity, are not already parametrized 
systems [23. For such models, there is no privileged time 
and no symmetry in general (cf. |29l|). But in asymptoti- 
cally flat cases, there is a privileged asymptotic time and 
an asymptotic symmetry. As it is shown in p^ . such 
symmetries do not determine their generators in the ex- 
tended phase space V uniquely but only up to addition 
of a linear combination of constraints. Despite that, they 
still define a unique dynamics of Dirac observables. 



III. POLARIZED CYLINDRICAL WAVES: 
SOLUTIONS AND ASYMPTOTIC BEHAVIOR 

A vacuum spacetime describing cylindrical gravita- 
tional waves with a fixed state of polarization (one degree 
of freedom per point) has two commuting, hypersurface- 
orthogonal spacelike Killing vectors d/d(p and d/dz; the 
Killing field d/dip is rotational and it keeps a timelike 
axis fixed; d/dz is translational; coordinates p> and z are 
invariantly defined up to a translation z l—^ z + a. The 
metric can be written in the form 

where T and R are invariantly defined, T up to a transla- 
tion T i—f T + a. In the above equation, ip = 'ip{T,R) and 
7 = 7(T, P). (To obtain Eq. (|3.1|) one uses a consequence 
of vacuum field equations — see, e.g., 0, [sof.) 

It is well-known that because of the translational sym- 
metry d/dz, the four-dimensional Einstein equations are 
equivalent to the three-dimensional Einstein equations 
with certain matter sources (see, e.g., |23| and [T^). 
In our case of cylindrical symmetry {d /dp is a further 
Killing field) the four-dimensional Einstein vacuum equa- 
tions the solutions of which give Einstein-Rosen waves 
are equivalent to Einstein equations in three dimensions 
with a zero-rest-mass scalar field ?/; as a source. It is, how- 
ever, more advantageous for the canonical formulation to 
work with the physical Klein-Gordon field = ^/VSG, 
G being the Newton constant. Hence, we formulate ev- 
erything with the help of the field (j). 

In three dimensions, the metric is given by (cf. and 

113) 

ds^ = gabdx''dx'' = e''{-dT^ + dR^) + R^dp'^ , (3.2) 
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and the Einstein field equations become 



u 7 


u 7 


1 8^ 


dB? 






927 


927 


1 97 












1 97 









= 8G 



9T 



8G 



d(j) def) 
dRdf 



The field equation for 

QT2 + 



1 



(3.3) 
(3.4) 
(3.5) 

(3.6) 



is the wave equations for the nonflat metric (|3.2|l as well 
as for the flat (Minkowski) metric obtained by putting 
7 = in (|3.2II . This crucial simplification implies that 
the scalar field is decoupled from the equations satis- 
fied by the metric. Eqs. (|3.3|I - H3.5|I reduce to two simple 
equations 



dT 



90 
OR 



of - ^^^dfdR 



(3.7) 
(3.8) 



the wave equation (|3.t)|) is their integrability condition. 
We can thus solve the axisymmetric — in three dimen- 
sions 'spherically' symmetric — wave equation (|3.6I) on 
Minkowski space and then solve Eqs. H3.6|l and l|3.8|) for 
the metric function 7(T, K) by quadratures. These well- 
known facts are of key importance in the canonical and 
quantum theory since all physical degrees of freedom are 
contained in the scalar field. 

We shall now briefiy review some of the results on the 
asymptotics obtained in and [sj. We shall extend 
the discussion by including both future and past null in- 
finities, and later also by employing a Fourier- type de- 
composition. 

The Cauchy data for the scalar field 0, given on the 
Cauchy surface topologically ^ suffice to determine the 
whole spacetime metric. For data which fall off appro- 
priately, the three-dimensional Lorentzian geometry is 
asymptotically flat both at spatial l2ll| and null infin- 
ity |23 | although in four dimensions the Einstein-Rosen 
spacetimes are not asymptotically flat (see for a de- 
tailed investigation of cylindrical waves at null infinity in 
four dimensions). 

By employing the "method of descent" from the Kirch- 
hoff formula in four dimensions one can find the represen- 
tation of the solution (/)(T, R) of the wave equation (|3.6() 
in three dimensions in terms of Cauchy data 0o — '/'(O, R) 
and 01 — 4>,t{0, R)- This has been used in to find the 
asymptotic behavior of the field and the whole metric 
(|3.2() at the future null infinity for the data of compact 
support (see Sec. II in ^^). By applying the same proce- 
dure one can analyze the solutions at the past null infin- 
ity. Introducing retarded and advanced time coordinates 



U = T- R , V = T + R 



(3.9) 



(notice that these are null coordinates for both fiat 
Minkowski metric and the curved metric (|3.2|) ). one ob- 
tains expansions in the powers of R~^/'^ along null hy- 
persurfaces U = constant and V = constant of the form 



0(^,i?) 
0(C/,i?) 



00 



9k{V) 
jlk+1/2 ' 



(3.10) 



;^/(^)-^E|g4- (3.n 

^ k=l 



The coefficients in the expansions are determined by the 
Cauchy data. By rewriting the Einstein field equations 
(|3.7|) and (|3.8|) in terms of U and R (respectively V and 
R), we obtain the asymptotic behaviour of the metric 
function 7 at 1+ in the form 



liU,R)^j. 



-8G 



dU 



dU 



+ 0{R-^) , (3.12) 



and similarly at . Here the constant 700 , which will 
play a key role in the following, is determined uniquely 
by the Cauchy data for (cf. Eq. (|3.7|l ) 



dT 



OR 



7oo = 7(0,00) = 4G / dRR 
Jo 

(3.13) 

The value of 700 represents the total energy of the 
scalar field computed by using the Minkowski metric. 
For any nontrivial data, 700 is positive. Hence, the metric 
at spatial infinity, given by 



ds' 



e^°° {-dT^ + dR^) -f R^dip^ 



(3.14) 



has a conical singularity because the distance of the cir- 
cles with radii R and R -\- dR is different by a factor e^^ 
from the difference of their circumferences divided by 2-k. 
It can be shown that as one approaches T'^ (i? ^ 00, 
U = constant), one finds (cf. Eq. (|3.12(l 'l 



7(f/,oo) =7oo-8G 



du ( 

du 



(3.15) 



and 7 to vanish at the timelike infinity Hence, 



7oo = 8G / dU 



(3.16) 



The conical singularity, present at spacelike infinity, is 
thus "radiated out", and the future timelike infinity i'^ 
becomes smooth. Eq. (|3.15|) plays the role of the well- 
known Bondi mass- loss formula, the function df /dU be- 
ing analogous to the Bondi news function (see also [s^ . 
Eq. (3.6), for an analysis in four dimensions). Clearly, 
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analogous formulae to (|3.15|) and (|3.16|) are valid for in- 
coming waves, with df/dU replaced by dg/dV: 

,iV,o.)^8G f dv 



and 



7oo = 8G/ dV{^ 



(3.17) 



Here we assume smooth past timelike infinity i~ and in- 
coming waves from the past null infinity 2^ with a null 
data g{V) bring in mass-energy which reveals itself as a 
conical singularity characterized by 7(V,cx)). At spatial 
infinity iq {V = oo, R = oo) this becomes just the con- 
stant 7oo given in Eq. H3.13|l . The fluxes of radiation, the 
analogues of the news function, as well as conical singu- 
larities are observable quantities at the past and future 
null infinities. Both are given by the asymptotic null 
data g{V) and f{U). The asymptotic null data will be 
important in the following. 

Starting from the representation of the solutions of the 
three-dimensional wave equation (|3.t)|) in terms of the 
Kirchhoff-type formula obtained by the "method of de- 
scent" from four dimensions one can, for the Cauchy data 
of compact support, obtain not only expansions (|3.1U|) 
and (|3.11|l . but also the explicit expression for the null 
data /([/) (respectively giV)) as the integral over the 
Cauchy data 0o and 0i. However, these integrals be- 
come simple only for retarded times U so large that the 
support of the data is completely in the interior of the 
pas t cone (similarly for the advanced times at see 
|23 |. Here we need the null data for all ?7's at and 
F's at 1^. 

To achieve this, we start from a Fourier-type decom- 
position. This, in three dimensions, means to write the 
solutions in terms of the Bessel functions of zero order 
provided that we require the solutions to be regular ev- 
erywhere, in particular at i? = (see, e.g., |33j). 

Thus, we start from the solutions of the form 



0(r,i?) 



1 



\/2 7o 



duo [A{u)jQ{ujR)t 



-iujT 



(3.18) 

as usual, we write just "c.c." instead of the second term, 
meaning the complex conjugate of the first one. Using the 
asymptotic expansion of the Bessel function at i? ^ oo 
(see, e.g., HJl), we obtain 



(/)(T,i?) 



1 



duj 



2\fRK Jo 



A{uj)e' 



-i{Tr /4)—iujU 



+ C.C. 



+ 0{R-^'^) , (3.19) 



where U and V are retarded and advanced time coor- 
dinates given by Eq. H3.9|l . Hence, the null data at the 



future and past null infinities read as follows: 



f{U) 



9{V) 



1 



duo r 



A{uj)e 



-i{n/4:)-iujU 



dhj r 



Jo 



A{uj)e 



(3.20) 



(3.21) 



It is easy to invert the last equations by writing, for ex- 
ample. 



giV) + 9{-y) = \ - I dw[A{uj)+ A* (oj)] COS ujV 



(3.22) 

g{V) - g{-V) = -i^ J duj [A{uj) - A*{lo)] sirnoV , 

(3.23) 

where A{uj) — {2Lu)^^/^A{uj)e^'^^^. Using Fourier cosine 
and sine (inverse) transforms to express A±A*, we find 



A{uj) = J -e-'^/^ / dV [.g(U)e^"^ + 9{-V)e- 



iujVl 



(3.24) 

A*{uj) being given by complex conjugation. Alterna- 
tively, we can write 



A{lu) = ^ ^e-"/4 J dV g{V)e'^^ 
Similarly, from Eq. (|3.2Q(I we get 



(3.25) 



A{u;) = ,/-e-/4 / dU [f{U)e'-^ + f{-U)e 



— iujUl 



or 



fiU)e 



(3.26) 



(3.27) 



Since according to Eq. H3.18|l the functions A{ijj) deter- 
mine the solutions 0(r, R) everywhere, Eqs. (|3.24(l - (|3.27|) 
imply that either the null data g{V) at I~ or f{U) at 1+ 
determine 4>{T, R) uniquely in the spacetime. 

The amplitudes A{lo) can be expressed also in terms 
of the Cauchy data 0o = '/'(0,i?) and 0i = (/)_T(0,i?) 
directly from Eq. H3.18(l . Using the Hankel transform 
(see, e.g., H^l): for two functions X{x) and Y{y), 

X{x)^ / dyY{y)^Jo{xy) (3.28) 
Jo 

is equivalent to 

Y{y)^ dxX{x)^Jo{xy) . (3.29) 
Jo 
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Expressing (f>i from Eq. H3.18|l we obtain 
1 

71, 



A{co) 



dR {lu(I)o - i(j)i)RJo{ujR) . (3.30) 



Hence, as expected, we need both 0o and (f>i to determine 
the solution of the wave equation (|3.6I) everywhere. For 
time-symmetric initial data, (f>i — 0, the amplitudes A(uj) 
become real. 

Although for the Cauchy data of compact support and 
even for more general data falling off sufficiently rapidly 
at spatial infinity we get ~ 1/^/R at null infinities as 
in Eqs. H3.10|l and H3.11|l (see [IJ), at spatial infinity, i.e. 
in the limit R oo, T fixed, the solutions fall off more 
rapidly: 



^^0{1/R), ^ 0(l/i?2) . 



(3.31) 



This will be needed in the canonical theory. To demon- 
strate the fall-off, employ again the asymptotic form of 
MljR) at i? c» in Eq. p!T5)l . 



const 



uj)e' 



and put ujR — lo' in the integral. Then we get 
const f°^' duj' 



: cos W 



-iuj'T/R 



which at large R and fixed T leads to ~ 1/ R. 

Finally, let us illustrate previous results b y a simple ex- 
ample. The Weber- Wheeler-Bonnor pulse [2J|, [231 rep- 
resents an exact, time-symmetric vacuum solution of the 
Einstein equations with cylindrical symmetry which sat- 
isfies all regularity conditions required above. The pulse 
comes in from the past null infinity, concentrates around 
the axis of symmetry (in three dimensions around the 
center i? = 0) at T = 0, and then reexpands to future 
null infinity. The real amplitude 



A{uj) = Ce 



(3.32) 



where C and a are constants, implies, by using Eq. H3.18|l . 
solution 



C 



[(a2+i?2_y2)2^4^2y2]l/2_ 



(a2 + i?2 - T2)2 + 4a2r2 



(3.33) 

regular everywhere. (Due to the factor 1 /-\/2 in Eq. (|3.18|l 
(j) here must be multiplied by \/2 to get Eq. (3.15) in [IJ.) 
At spatial infinity, R oo, T fixed, we see that 



evil 



0(l/i?2 



(3.34) 



in accordance with Eq. H3.31|l . At the past null infinity 

{R~^ oo, V = T + R fixed), we find 



of 



2U/2 



V2 



1/2 



(3.35) 



(At future null infinity, i?— >oo,l7 = T — R fixed, 
the same expression, with V replaced by U follows.) A 
simple calculation, starting from the formula H3.21(l for 
the profile g{V) and using A{uj) from Eq. (|3.32l) . yields 
exactly the factor at 1/^/R in Eq. H3.35|) (for integrals 
dx e~'^^{l/^/x) cosbx and dx e~'^^{l/ y/x) sinbx 
needed in the calculation, see e.g., [s^, formulae 3.944, 

13 and 14). 

With (j) given by Eq. 1)3. 35|) one can find the explicit 
expression for function 7 by solving Eqs. H3.7|l and (|3.8|l . 
It reads Hi, 113 as follows: 



^ \a2 [(^2 + ^2 „ y2)2 + 4a22^2]2 



1 



R' -a^-T^ 



(3.36) 



a2 [(a2^^2_y2)2_^4„2y2]l/2 

The conicity at spatial infinity is thus given by 

7oo = 8G(C/a)2 . (3.37) 

With this explicit solution one can verify directly rela- 
tions H3.15|l - I|3.17|l for the conicity as it is radiated "in" 
and "out" , similarly to the Bondi mass in four dimen- 
sions. 

IV. HAMILTONIAN FORMULATION 

Let us essentially repeat the canonical treatment of 
Ashtekar and Pierri given in 20] but with a slight modi- 
fication in order to establish the analogy with the model 
of Sec. 2. Let us start by considering the volume part 
of the canonical action derived from the Einstein-Hilbert 
action by assuming cylindrical symmetry in |20| : 



1 



S 

where 



8G.I 7 



dt I dr{p^j^t+PRRt+Pi^'4',t-NC-N'^ Cr) , 

(4.1) 



C= e-''/2 l^2R,rr--f,rR,r-p-yPR. + R~^pl/'2 + R'lp^r/'2) 



and 

Cr = e""^ ( - 2py^r + Py l,r + PR. R,r + Pil) V'.r) 

are the constraint functions, 7, R, and ip are defined by 
Eq. H3.1|l . p^,, PR and are the conjugate momenta, 
while N and N"^ are Lagrange multipliers — the so-called 
lapse and shift functions. One should add to this action 
the boundary energy term 



(4.2) 



and specify the fall-off of N accordi ng t o Ivair^oo N = 1 
so that the action be differentiable [2^. However, this 
term is not invariant under reparametrizations of the la- 
bel time t as there is no temporal density present in the 



9 



integrand in (|4.2() . Addition of the bare term (|4.2(l to the 
action (|4.1I) would imply a privileged choice of asymp- 
totic time. The total action would then not be in the 
constraint-Hamiltonian form but rather in an already re- 
duced form at spatial infinity. 

In order to recover the full constraint-Hamiltonian 
framework of our model in Sec. 2, we need to justify 
the inclusion of a temporal density in (|4.2|l . This can 
be done following the general approach by Beig and O' 
Murchadha "2^. First, one considers fall-off conditions 
for the configuration space data at r ^ oo. These have 
been specified in 20"|. The configuration space fields ap- 
proach infinity according to 



7(t,0->7oo(i) + 0(l/r) , 
R{t,r)^r{l + 0{l/r)) , 
^{t,r)^0{l/r) , 



(4.3) 



where r 0(l/r), O^i/r^), etc. admit limits at r ^ oo. 
(These limits generally depend on the time t.) One then 
requires that the action of the Liouville form 



dr 



{p^ d'y + pRdR+ dip^ 



on vector fields of the form 

should be finite. The resulting integral 

dr {p^ (57 +prSR + p^ Sip 



is finite if the momenta satisfy the following fall-off con- 
ditions at r ^ oo: 



p^{t,r) ^ 0{l/r^) , 
PR{t,r)^0{l/r^) , 
p4t,r) ^ 0{l/r) . 



(4.4) 



Next, concerning the behaviour of the lapse and shift, 
there are several aspects that ought to be kept in mind. 
First, it is the finiteness and differentiability of the Hamil- 
tonian part 



HilN, N'' 



dr {NC -f iVa) 



(4.5) 



of action glj (cf. g^). Second, Hi[iV,iV] has to gen- 
erate a transformation within the phase s pac e defined by 
the boundary conditions (|4.3|) and (|4.4|l 23]. Finally, if 
we are going to have a full analogy to action (|2.14|l of 
Sec. 2, we have to p arametrize the model also at infinity, 
as it is done in "24] for a spherically symmetric model. 

The constraints functional (|4.5(l remains finite even if 
N{r) and N^{r) approach arbitrary temporal densities 
at r — !■ cx); namely, 



7V(<,r)- 
N^{t,r) 



-0(l/r) 
-0(l/r) 



(4.6) 



where Noo{t) has to be non- negative. The condition that 
the lapse and shift should approach temporal densities 
at r — > cx) is the minimum requirement that is com- 
patible with the invariance of the action H4.1fl under 
reparametrizations of t. 

Conditions (|4.3|I - H4.6|1 now imply that the action is not 
differentiable. The problem comes from the variation of 
7(t, r) leading to the boundary term 



SS 



= J di(iVooe-^~/2 57oo) (4.7) 

at spatial infinity. In order to have a consistent canonical 
theory, one needs to add to the action the boundary term 

whose variation with respect to "f{t, r) cancels the bound- 
ary term in 14. 7|) . The differentiable action is therefore 



= ^ y" y" dr{p^ 7,t -I- PR R^t + V',* 

-NC-N''Cr) (4.8) 

The boundary term in 14.8|l has now been modified by the 
addition of a constant in order that it coincides with the 
asymptotic energy derived from first principles in [20l| . 
There is no boundary term involving the shift in spite of 
the fact that its asymptotic value may be non-zero. The 
corresponding constraint functional generates an "even 
supertranslation" (in the language of ^3|) and is differ- 
entiable without any correction, similarly to the case of 
four-dimensional general relativity, cf . 23] . 

When varying the action H4.8|l with respect to the lapse 
N(t,r), it is important to keep the ends of its variation 
fixed. Indeed, if Noa{t) is varied in (|4.8I) then one gets 
an unwanted field equation implying that the asymptotic 
energy vanishes. It follows that the action H4.8|l is not yet 
in true constraint-Hamiltonian form. Following Kuchaf 
j24j . this can be improved by the "parametrization at 
infinity" : iVoo {t) should be replaced by a differentiated 
asymptotic time dTao/dt = Too.t{t)- The asymptotic time 
is determined by the asymptotic metric: it must hold 
Noo = 1 if the parameter t is chosen to be Too ■ The time 
Too{t) can be varied in the ensuing action. Its variation 
leads to a redundant equation amounting to the conserva- 
tion of the asymptotic energy. One should next introduce 
the momentum P^xd and add the associated constraint 
(which is linear in Poo) to the action by a new Lagrange 
multiplier Noo{t). 

In this way the action is brought into the true 
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constraint-Hamiltonian form 



S — / dt ( Poc Too,t 



+ y rfi y" dr (py 7^t + PR Rt + P-^ 
- JdtN^ (Poo + ^ (1 - e-^-/2)) 

JdtJ^ dr [n C + N"- Cr) . (4.9) 

The multipliers N and N"^ obey the asymptotic condi- 
tion (|4.6|l . The action (|4.9(l is the analogue of the action 
(|2.14|) for the Newtonian particle. One can verify that the 
field equations derived from the variations of 14. 9() coin- 
cide with those of Sec. 3, preserve the fall-off conditions 
(|4.3|) - (|4.()|l and imply the conservation of the asymptotic 
energy. 

Action (|4.9() is our starting point for the canonical the- 
ory. Although it corresponds to action (|2.14|) of Sec. 2, 
observe that there is no a priori analogue of action (|2.3(l 
of Sec. 2. We have to begin with the extended phase 
space V with coordinates ^(r), p-y{r), R{r), pnir), ip{r), 
p^{r), Too, Poo and the symplectic form fl-p, 

1 1 f°° 

+ dpR{r) A dR{r) + dp^{r) A d^{r)\ . 

The constraint surface C is defined by TilN^N"^'] = 
for all N{r) and N^ir) satisfying the fall-off conditions, 
where 



n[N, = TV^ ( + — (1 - e-Toc/2; 



— y dr{NC^N^Cr) 



(4.10) 



The orbits are defined by the canonical action of the con- 
straint functional 7i[A^, A^*"]. The canonical transforma- 
tions generated by H4.1()|l are considered as a gauge trans- 
formation. Within the gauge group, there is no distinc- 
tion between the "symmetry" and the "proper gauge" as, 
e.g., in [23, m and . The functional (jCTTjl gener- 
ates only reparametrizations both "inside" the spacetime 
and at infinity. Symmetries are now generated by differ- 
ent functions. Indeed, the functional H4.10I) has vanishing 
Poisson brackets with Pqo for any N{r) and N^{r) sat- 
isfying the conditions (|4.6|l because the variable 700 is 
asymptotic value of the canonical coordinate 7 so that 
{-PoojToo} = 0. Hence, it is the function Poo that gen- 
erates the symmetry. In general, we conjecture that one 
can introduce privileged coordinates at infinity and that 
asymptotic symmetries are generated by their conjugate 
momenta or suitable combinations of the momenta and 
the coordinates (like, e.g., boosts). 

The variable Too to which Poo is conjugate is a kind 
of a "privileged time" but the surface Too = constant is 
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FIG. 1: Important surfaces in the constraint manifold C. 
The intersection of Too = constant with any orbit is infinite- 
dimensional. The gauge condition surface Q intersects each 
orbit in a (one-dimensional) dynamical trajectory of the re- 
duced theory. The points common to Q and each surface 
Too ~ constant is a (infinite-dimensional) transversal surface 



neither a transversal surface in the phase space, nor a 
Cauchy surface in each solution spacetime. Indeed, the 
function Too— constant has vanishing Poisson brackets 
with HlN, N"^] for all N{r) and N^{r) whose asymptotic 
values vanish; hence, the duly generalized regularity con- 
dition 12.19|l is not satisfied. It follows that an infinite- 
dimensional submanifold of each orbit lies in the surface 
Too — constant (Fig. 1). This is connected to the fact 
that the condition Too = constant defines only a particu- 
lar section of infinity in each cylindrical wave spacetime 
but not a Cauchy surface of the whole spacetime; there 
is a relation between Cauchy and transversal surfaces, cf. 

The reduction by gauge condition, analogous to that 
described in Sec. 2.1, starts by a choice of a one- 
dimensional family of transversal surfaces. Let us denote 
the manifold formed by all chosen transversal surfaces in 
C by G- In Sec. 2, a privileged choice of gauge has been 
possible: Q has been the family of surfaces T — t, t E M., 
where T is the privileged time. The nearest to this we 
can come is to choose the transversal surfaces in Q to be 
the intersections of G and Too = constant (Fig. 1). There 
are, of course, many choices of Q. One example of such a 
choice is carried out in . Let us describe an analogous 
choice for our action H4.9|l . 

Following Ashtekar and Pierri, one may fix the part of 
the gauge associated with the constraints C(r) =0 and 
Cr{r) = by imposing the gauge- fixing conditions 



R{r) = r , p-y(r) — . 



(4.11) 



These are the defining equations for Q. Viewed as con- 
straints, these conditions form together with the con- 
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straints C(r) =0 and Cr{r) = a second-class system. where 
The remaining constraint Poo + 47; (-'^ ~ e~'''°°/^) — can 



be taken care for by the gauge-fixing condition 



constant 



(4.12) 



The surface Tt^ in C defined by H4. 1111 - 14. 121) selects an 
initial datum from each gauge orbit in C. The gauge- 
condition surface Q is swept by all Tt^ . 

In order to confirm that this reduction is admissible, let 
us add l|4.1HI and (|4.12() to the action H4.9|l by Lagrange 
multipliers M, Af, and find out if the ensuing action 
determines unique values for N, TV''. One obtains the 
action 



1 



S — / dt ( Poo Poo,t , 



^Poc Poo,i^ 



1 



4G 



l-e 



dtMoo{Too-t) 
l,Jdtfdr{NC + N^a) 



dt 



driM (R 



(4.13) 



where the set of conditions (|4.12(l is implemented by the 
expression J dt M^c (Poo — t) . Indeed, it is not difficult 
to check that all redundant variables in (|4.13() can be 
expressed uniquely in terms of the canonical pair {ip,p,p) 
by solving the set of equations derived from the variation 
of these variables in (|4.13l) . This confirms that the gauge- 
fixing conditions 14.11|l arc regular. In particular, the 
unique reduced expressions for the multipliers N, and 
Noo are 



^(Poo,P) 

= exp 



dr 



iV^(Too,P) = 



(4.14) 



those for the canonical pairs (7,^7), {R,pr) and 
(Poo,PToo) read 



7(^00, P) 

p^{Too,R) 

P(roo,r) 

Pr{Toc,R) 
Too{t) 



dr 



{r-'pl{T^,r)+ri,^,{T^,r)') 



2 Jo 

, 

r , 

-P-0(Poo,P) V',-R(T'oo,P) , 
t , 



Poo (Poo ) — 



1 - exp ( - ^700 



7o 



dp(p-14(Too,P)+PV''fl(Too,P)) , 



(4.15) 

and unique expressions also follow for the multipliers M, 
AP, Moo- The uniqueness of these expressions partially 
relies on the conditions imposed on the canonical fields 
at r = (see, e.g., |2(|) which force 7(Poo,0) to vanish 
for aU Poo. 

The reduced action for the remaining canonical pair 
(■0(Poo, P),P0(Poo, P)) on Q C C, parametrized by the 
values of the asymptotic time, is therefore 



5=^y"dPoo^ dR[p4T^,R)i^,T^{To,,R) 



AG 



dT^ 



(l„e-^-(^-)/2) , (4.16) 



where 7oo(7oo) is expressed as a functional of ^(Poo,P) 
and p^{Too,R)) in H4.15|l . The action (|4.1t)|) is analo- 
gous to the reduced action H2.3|l . The phase space Fi 
is described by coordinates i^iR) and p^{R), while the 
symplectic form is 



dRdp^{R) Ad^{R) 



The action 14.16|l is precisely the reduced action of 
Ashtekar and Pierri. In particular, the Ashtekar-Picrri 
time t in Eq. (19) of |23| corresponds to the time Poo 
here. 

Geometrically, ipir), p^{r) and Poo are coordinates on 
the gauge-condition surface Q. The surfaces defined in Q 
by the equation Poo = constant are transversal surfaces in 
the phase space V . However, they also determine a fam- 
ily of Cauchy surfaces of constant Ashtekar-Pierri time in 
each solution spacetime (cf. H^])- In this sense, the part 
(|4.11|) of the gauge condition determines a particular ex- 
tension of the points at infinity defined by Poo = constant 
to whole Cauchy surfaces in the spacetimes. However, 
different choices of G lead to different Cauchy surface ex- 
tensions of these points at infinity. Hence, two different 
choices of Q entail two different choices of time so that the 
transformation between the times has again the character 
of Eq. H2.20|l even if the part (|4.12() of gauge conditions 
remains always the same — only the asymptotic values of 
these times have then to coincide. As noted at the end 
of Sec. 2.1, it is the transformation (|2.2U|) between re- 
spective times which causes difficulties in constructing a 
unique plausible quantum theory. 

Considering the privileged symmetry generated by Poo , 
we can see that it remains a symmetry of the reduced 
theory. It acts in Q as follows. 



[<Pir),p^^{r),Too^ (^i/'(r),p^(r). 



P« 



(4.17) 



12 



while the original action of Poc in V is 



(7('') , Pj (r) ,Rir), PR (r) ,ip{r), (r) , T^, 

^ {l{r),Pfir),R{r),PBir),Tp{r),p^{r),Too +■ 

(4.18) 

The map (|4.18|l is tangential to Q and the map (|4.17|) is 
just the restriction of (|4.18|l to G- This follows from the 
fact that the constraints as well as relations H4.11|l that 
define G are independent of Too- 

The dynamics defined by action H4.16|l determines a 
foliation of G by one-dimensional dynamical trajectories 
represented by two functions of two variables iplRjToo) 
and p,p{R, Too ) ■ These are identical with the intersections 
of G with the orbits. In this way, we obtain a bijection 
between integrals of motion of the reduced theory and 
Dirac observables. On one hand, any Dirac observable is 
constant along each orbit. Hence, it must also be con- 
stant along each dynamical trajectory of action H4.16|l . 
On the other, any function on G that is constant along 
each dynamical trajectory defines a unique extension to 
C that is constant along each orbit. 

This relation between the Dirac observables of the ex- 
tended system and the integrals of motion of the reduced 
theory, together with the compatibility of the symmetry 
groups generated by Poo in the extended and reduced 
theories, justify the approach of Sees. 5 and 6, where 
we shall construct the gauge-invariant dynamics starting 
from the gauge-dependent action H4.16|l . 



V. PHYSICAL PHASE SPACE Ta 

In this section, we choose a complete set of Dirac ob- 
servables, find their Poisson algebra and calculate their 
Poisson brackets with the symmetry generator Poo . This 
task is simplified if we start from Ashtekar-Pierri reduced 
action (|4.16|) instead of the original parametrized action 
(|4.9() . According to what has been shown in the previous 
sections, the result is independent of the gauge chosen to 
reduce the action H4.9|) . 

The reduced action (|4.16l) can be rewritten in terms 
of the rescaled field — V'/VSG introduced in Sec. 3 as 
follows: 



S = J dtdR (7r00 - H) , 



where 



1^=^G I dR[ -n'^ + Rct>'- 



enters the Hamiltonian 



(5.1) 



(5.2) 



For simplicity, the notation for our time Too has been 
changed to the Ashtekar-Pierri notation t. The Hamilto- 
nian depends on t only through tt^ and so that H and 
7oo are constants of motion, 



7o 



. 



(5.3) 



The canonical equations that follow from the action 
are 



7o 



fit — f=^2^ 



■{Rcf^y , 



(5.4) 
(5.5) 



Eqs. (|5.5|) and H5.3() imply 



so that 



e 2 



9^ 



R 



(5.6) 



If we use the relation between the Einstein-Rosen time T 
and the Ashtekar-Pierri time t (see poj'). 



T{t) = e~^^-^t 



(5.7) 



then Eq. H5.(j|l becomes the wave equation (|3.t)|) . The 
general solution to Eq. H3.6|l is given by Eq. H3.18|l . which 
can be written in terms of time t as 



cf){t, R) 



1 



and Eq. (|5.5|l yields 
R 



duj 



dhJ 



A{Lo)Jo{uoR)e~'^^'^'^ 



(5.8) 



-it^A(tj)Jo(wi?)e-*"^(*) +. 



(5.9) 

Eqs. (|5.8() and (|5.9() describe the general solution to the 
canonical equations H5.4|l and (|5.5I) in terms of the set of 
constants Aiuj). Hence, the parameters A{lo) can serve 
as coordinates on the physical phase space 

The physical phase space is a symplectic manifold. Its 
full structure can be obtained if we find a transversal 
surface. As has been explained in Sec. 2, any transversal 
surface, together with the symplectic form that results 
from pulling back the symplectic form from the extended 
phase space to the transversal surface, form the struc- 
ture that is isomorphic to the physical phase space. In 
our case, the initial data (p^ and tt^q of the canonical coor- 
dinates (j) and TT^ at the Cauchy surface t = Q determine 
a unique solution (|5.8|l and (|5.9|l so that they can also 
be considered as coordinates on the physical phase space 
Fi. Moreover, the surface Tq defined by the Ashtekar and 
Pierri gauge H4.11|l together with the condition Too — 
is a transversal surface. Hence, the symplectic form on 
the physical phase space with respect to the coordinates 
00 and TT^o is 



dRdn^oiR) A d0o(i?) 



(5.10) 
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because this is the pull back of ilp to Tq by the injection 
map of Tq into V] the manifold Tq with this symplectic 
form is isomorphic to the physical phase space 

The relations between the parameters A{lli) and 0Oi 
TT^o can be obtained from Eqs. H5.8|l and H5.9|l : 



MR) = ^J duJo{LuR)[A{Lu)+A*{Lo)] , (5.11) 



and 



iR 



Tr^o{R) = -^ I dLULuJoiLoR)[A{Lu)-A*{oj)], (5.12) 
V 2 Jo 

while the inverse transformation is analogous to Eq. 



A{lu) = ^ / dRJo{LoR)[ojRMR) - i^^R)] ■ 
Jo 

(5.13) 

A further set of parameters to determine points of the 
physical phase space are the T~ null data g{V) or X+ 
null data f{U). The transformations between A{u!) and 
g{V) is given by Eqs. H3.21|l and H3.25|l . those between 
A(w) and f{V) by Eqs. (lOHl and iTOTIl . 

The quantity 700 is a function on the physical phase 
space given, in terms of the four different coordinate sys- 
tems, by Eqs_ETj, (|XTB|) a nd (|XT7|) . Eq. (|5TJ, into 
which Eqs. H5.11|l and (|5.12(l are substituted, yields af- 
ter some simple transformations the expression for in 
terms of A{uj): 



j^^SG I dujLuA*{uj)A{Lu) 
la 



(5.14) 



We can also express the symplectic form ()5.1()(l in terms 
of A{uj). If Eqs. H5.11|l and (|5.12(l are substituted into Eq. 
(|5.1U|) . we obtain 



n 



-- / duj j dcj' dRLj'RJo{ujR)Jo{uj'R) 
2 Jo Jo Jo 

X [- dA{uj) A dA{uj') + dA{uj) A dA*{uj') 
- dA* [uj) A dA{uj') + dA* (uj) A dA* (w')] 
The formulae (|3.28|) and H3.29|l imply, however, that 



/ dRRJo{ojR)MLu'R) = -^=SiLu - uj') . (5. 
Jo Vluuj' 
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Hence, using the antisymmetry of the wedge product, we 
obtain finally 



n2^i dujdA*{uj) AdA{uj) . (5.16) 
Jo 

Let us also express the symplectic form of the physical 
phase space in terms of the asymptotic null data g{V) 
and /([/). Eqs. fCTH^ and (TT^ give 



Since the wedge product is antisymmetric in V and V', 
only the antisymmetric part of the integral over uj con- 
tributes to the result. However, 



— / dujuj 
2vr7o 



1 

2^ 



1 d 
2^dV 



Hence, 

/CO 
dV dg'{V) hdg{V) . (5.17) 
-00 

By analogous calculation, Eqs. (|5.16|) and (|3.27|1 yield 

/oo 
dU df'{U) hdf{U) . (5.18) 
-00 

Finally, let us calculate the transformation between 
f{U) and g{V) if f{U) is defined by the solution de- 
termined by g{V). Such a transformation is, therefore, 
entailed by Eqs. and ^Z^: 



dVg{V) / duj 

a Jo 



^^iuj{U-V)_^^-tu{U-V) 



The distribution D{U — V) defined by the integral over 
UJ can be approximated by a convergent series of distri- 
butions - V) (see l37j), 

\iiaD^{U -V) = D{U~V) , 

e^O 



where e > and 



D,iU-V) 

Jo 

= -2 



However, 



duj 



u-v 



lim 

£-►0 



dV 



U-V 



{U - Vf + £2 
= -2P 



9{V) 



dV 



9{V) 

u~v ' 



where P denotes the principal value. Hence, 

9iV) 



fiU) 



1 

-P 

TT 



dV' 



u-v 



(5.19) 



VI. REPRESENTATION OF SYMMETRIES IN 
THE PHYSICAL PHASE SPACE 



i / dujdA*{uj) A dA{uj) 



I 

TT 



-00 J —00 



dV I dV dg{V)hdg{V) j duj uje"^^^-^'> 



There are two interesting symmetries to be represented 
in the physical phase space. The first is the infinitesimal 
time translation, and the second is the map 
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defined by ?7 = V" in terms of coordinates U at X"*" and V 
at T~ (an analogous symmetry transformation has been 
studied in 25]). 

Tlie push-forward action of the infinitesimal transla- 
tion t^t + dtow the solution fields (l){t, R) and 7r0(t, R) 
is given by 

(l){t,R) ^ (j}{t-6t,R) , (6.1) 

TT4,{t,R)) ^ TT^{t-St,R) . (6.2) 

Substituting t — St for t into Eqs. (|5.8|l and H5.9|l and com- 
paring the results with these equations leads to A{lu) 
A(lo) + SA{u)), where 

6A{uj) = iLue-^^°°A{Lj)St . (6.3) 

The same result can be obtained, if we put 4>o{R) + 
d(j)o{R) and tt^o{R) + STr^o{R) into Eq. I|5.13(l and calcu- 
late the corresponding SA{uj). We must utilize the fact 
that 

express the time derivatives with the help of the equa- 
tions of motion H5.4|l and H5.5|l , transfer the r-derivatives 
from (j)o{R) to Jo{lljR), and use the Bessel equation that 
is satisfied by Jo{ujR), 

-^{R4{ujR))' = uj^UujR) . 

It follows that the action of the infinitesimal time 
translation is canonically generated by the function —H 
defined by Eq. (|5.2|l . with given by Eq. (|5.1(l . where 
(j) and 7r0 are replaced by <j)t^ and tt^q. (Indeed, the mo- 
mentum conjugate to t is Poo = .) We thus have 

STT^oiR)^{7^MR)^-H} , 

SAiio) ^ {Aiu;)-H} , 

and obtain analogously 

Sf{U) = {/([/), -i/} , (6.4) 
SgiV) = {giV),-H} . (6.5) 

Let us express Sf{U) and Sg{V) explicitly from the ac- 
tion of translations. Since the whole solution is shifted 
along the background manifold defined by the coordi- 
nates t and R hy 1 1-^ t + dt, R 1-^ R, we have, regarding 
Eq. CT . 

U ^T- R^ e"^^~ {t + 6t)- R=U + e-^''°-St , 
and similarly for V: 

V + e^i'''-^St . 

Hence, 

Sf{U) = -f'{U)e--^^-5t , (6.6) 
5g{V) = -g'{V)e~^'<^5t . (6.7) 



The same relations result from the Poisson brackets (|6.4(l 
and inSIl, if Eqs. JX^ . fOIjl and (jTTCI) 

are used; notice that Eq. (|5.16|l implies 

{A{Lu),Aico')}^0 , {A*{lu),A*{lu')}^0, (6.8) 

and 

{A{u;),A*{lu')} = -iS{cj-uj') . (6.9) 

As has been explained in Sec. 2 (see also 0), the dy- 
namics of the Dirac observables is defined by the compar- 
ison of the equations of motion with the action of symme- 
try. Since the equations of motion for the Dirac observ- 
ables are trivial (the observables remain constant), the 
symmetry action alone gives the total dynamical change. 

The second symmetry a : 2^ t—f is a purely asymp- 
totic one, similarly to Too ^ Too + t. Its action on solu- 
tions can be found in an analogous way via the Cauchy 
data for solutions at I's. We consider Cauchy null datum 
gi{V) at T~ as defining solution (/'i(t, R). Then we push 
forward the field gi{V) at to X+ by cr*, which results 
in Cauchy datum f2{U) at X+. The Cauchy datum f2{U) 
determines another solution 4'2it, R), and we define it as 
the image of 0i(i, i?) by a. The corresponding map in 
r2 can be calculated by using coordinates f{U) in 
Solution (j)2{t,R) has coordinate J^iU), and let (pi{t,R) 
has coordinate fi{U). Then the point f2{U) in r2 is the 
image of the point fi(U) of r2 by map a. The dynamics 
defined by a is the inverse map because it compares the 
evolution by the wave equation (which is trivial because 
Dirac observables remain constant) with the map by a 

(cf. Hi). 

The push forward map of fields at I to those at 
1+ acts as follows: 

a.g{V) = f{U) , 

where f{U) — g{U). It follows immediately that the 
dynamical evolution defined by the "zero motion" a is 
represented by transformation (|5.19|) . 

We can also introduce Fourier amplitudes a{uj) and 
b{uj) of the asymptotic data by 

a{Lu) = yl(cj)e*^/'* , b{Lu) = ^(t^)e-"/4 , 
so that Eqs. I|3.20(l and H3.21|l become 



f{U) = 
9{V) - 




The push forward of the amplitudes is clearly given by 

cr*a(cj) = b{uj) , 

where b{uj) — aluj). Canonical representation of a is, 
therefore: 

b{uj) = a{ij)e-'''/^ = ~ia{uj) . (6.10) 

This map (or H5.19|l 'l becomes the S'-matrix of the one- 
particle sector in the quantum theory of the model. 
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VII. QUANTUM THEORY 

It is easy to construct the Hilbert space, the operators 
representing the Dirac observables, the Hamiltonian, and 
to define the scattering matrix in the standard way of 
quantization of Hnear field theories (see, e.g., or [s^V 
A sketch thereof wih be described in this section. 

Let us start from the Poisson brackets (|6.8I) and (|6.9() 
for the observables A{uj). Roughly, in the canonical quan- 
tization, Poisson brackets are replaced by commutators 
multiplied by i (the units are chosen so that the Planck 
constant is 1). Then, we have 

[A{ij'),A{u;)]=0 , [A{ij'),A\u;)]=S{uj' -uj) . (7.1) 

These are commutators of the annihilation and creation 
operators of a quantum field theory for a continuous spec- 
trum. They form our starting point. 

For many constructions it is favorable to use a smeared 
version of the operators. We choose any complete or- 
thonormal basis of (complex) functions X„{lu), where 
Lu € (0,oo). This means that any complex function/ 
can be decomposed. 



/(u;) = ^/„X„(c^) , 



where fn are complex coefficients, and that 

/ duj X:^{uj)X,n{Lu) ^ 5nm . (7.2) 

Jo 

Defining 

i„= / duj X:^{uj)A{uj) , (7.3) 
Jo 

we obtain 

i(tu) = ^X„(w)i„ , (7.4) 

n 

and 

Then we can define the vacuum state |0) by 



i„|0) =0 Vti , (0|0) 1 
which also implies that 

i(w)|0) =0 . 



(7.6) 



(7.7) 



The elements of a complete basis in the Hilbert space 
are obtained by application of any number of creation 
operators A}^ to |0); if the total number of the creation 
operators is N, then the state is an 7V-graviton state. 
The scalar product is defined by scalar products of the 
basis elements, which, in turn, are determined by the 



commutation rules 1)7. 5|l and the conditions 1)7. 6(1 . For 

example, 

(iLlo),iJ.|o)) = (o|i™it|o) 

= {0\Al^An + 6nm\0) = S„rn ■ 

The Hilbert space defined in this way is often called the 
Fock space and we denote it by J-. 

Those Dirac observables defined in Sec. 5 that are lin- 
ear in the variables A{(jj) and A'^{oj) can be associated 
with operators on that are linear combinations of the 
operators A{uj) and A''{uj) with the same coefficients. 
This definition preserves the relation between Poisson 
brackets and commutators. For example, we define 



f{U) 



1 



dw 



A{uj)e 



— i{-!T / i)—iu)U 



+ it(c^)e*(-/4)+-t^ 



The matrix elements of f{U) with respect to the Fock 
basis are easily calculated by using the decomposition 
(|7.4|) . In such a way, we have a Hilbert space and the 
operators that correspond to the basic quantities. 

In order to construct the Hamiltonian, we start from 
Eqs. (|5.2|) and (|5.14|) . We define the quadratic operator 
7oo by the normal factor ordering: 



7oo = 8G / dujujA\ijj)A{uj) = 
Jo 



nm 



^nniA^nAra ; 



where 



/•oo 

8G / dwwX;(w)X„(w) 



Then 7oo|0) = 0. The operator 700 is self-adjoint on J^] 
it has a continuous spectrum. Its (generalized) eigenvec- 
tors form a (5-normalized basis of T , elements of which 
are obtained from the vacuum by application of any num- 
ber of the creation operators ^^(a;) (and a normalization 
factor). For example, 

7oo (it(c^)lO)) =8G^ {A\uj)\Q) 

Then, any function of 700 can be defined by the spectral 
theorem (see, e.g., ^as the same eigenvectors, 

and its eigenvalues are the values that the function has 
on the corresponding eigenvalues of 700 . In this way, the 
Hamilton operator 



H 



1 

4G 



1 — exp 



1 



is well-defined. For example, 
1 
4G 



1 — exp 



77o 



X 







because 700 has the eigenvalue zero on |0). 

Finally, we can define the scattering matrix S. In or- 
der to do that, we have to determine what are the in- 
and out-states. It seems natural to take the states that 
result from applying any number of the operators a) {uj) 
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to |0) corresponding to the observables a(w) of Sec. 6 
as the in-states. Shmlarly, the out-states can be defined 
by b{u!). From Eq. H6.1()|l . we have a simple Bogolyubov 
transformation between a(w) and b{u}): 

a{uj) — ib{uj) . 

The construction of the scattering matrix that imple- 
ments a give n Bogolyubov transformation is described 
in or |38l | . We shall skip it because it lies outside the 
scope of this paper. 

Another interesting question is, what is the relation 
between the Hamiltonian H and the scattering operator 
S. There are methods of calculating S from H: one has 
to take some limits within the Euclidean regime (see, e.g.. 



|13|'). However, an application, or even an applicability, 
of these methods to our case also lies outside the scope 
of this work. 
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